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Abstract 

In this paper we apply the techniques of symbolic dynamics to the analysis of a 
labor market which shows large volatility in employment flows. In a recent pa- 
per, Bhattacharya and Bunzel [1] have found that the discrete time version of the 
Pissarides-Mortensen matching model can easily lead to chaotic dynamics under 
standard sets of parameter values. To conclude about the existence of chaotic dy- 
namics in the numerical examples presented in the paper, the Li-Yorke theorem or 
the Mitra sufficient condition were applied which seems questionable because they 
may lead to misleading conclusions. Moreover, in a more recent version of the pa- 
per, Bhattacharya and Bunzel [2] present new results in which chaos is completely 
removed from the dynamics of the model. Our paper explores the matching model 
so interestingly developed by the authors with the following objectives in mind: 
(i) to show that chaotic dynamics may still be present in the model for standard 
parameter values; (ii) to clarify some open questions raised by the authors in [1], 
by providing a rigorous proof of the existence of chaotic dynamics in the model 
through the computation of topological entropy in a symbolic dynamics setting. 

Key words: Symbolic Dynamics, Periodic Orbits, Chaos Conditions, Backward 
Dynamics, Matching and Unemployment. 



Introduction 



In a recent paper, Levin [11] discusses the set of results presented over the last 
decade by various prominent physicists which led to the conclusion that black 
holes seem to be susceptible to chaos. Levin argues that the most realistic 
description available of a spinning pair of black holes is chaotic motion, and 
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goes on to complain that in physics and cosmology "chaos has not received the 
attention it deserves in part because the systems studied have been highly ide- 
alized". In contrast, in economics we have the interesting fact that even some 
of the most simple and highly idealized models describing modern economies 
can easily lead to chaotic dynamics. 1 

In this paper we apply the techniques of symbolic dynamics to the analysis of 
a labor market which shows in almost developed economies large volatility in 
employment flows. The possibility that chaotic dynamics may arise in modern 
labor markets had been totally strange to economics until recently, at least as 
far as we are aware of. In an interesting paper, Bhattacharya and Bunzel [1] 
have found that the discrete time version of the Pissarides-Mortensen match- 
ing model, as formulated in Ljungqvist and Sargent [12], can easily lead to 
chaotic dynamics under standard sets of parameter values. However, in order 
to conclude about the existence of chaotic dynamics in the three numerical 
examples presented in the paper, the authors apply the Li-Yorke theorem or 
the Mitra sufficient condition which should not be generally applied to all spe- 
cific simulations because they may lead to misleading conclusions. Moreover, 
in a more recent version of the paper, Bhattacharya and Bunzel [2] present 
new results in which chaos is completely removed from the dynamics of the 
model. This paper explores the matching model so interestingly developed by 
Bhattacharya and Bunzel with the following objectives in mind: (i) to show 
that chaotic dynamics may still be present in the model for standard param- 
eter values, for high values of the measure of labor tightness (that is the ratio 
of vacancies to the number of workers looking for jobs) which can occur in 
economic booms; (ii) to clarify some open questions raised by the authors 
in their first paper, by providing a rigorous proof of the existence of chaotic 
dynamics in the model through the computation of topological entropy in a 
symbolic dynamics setting. 

Therefore, if one is studying whether there are chaotic dynamics or not under 
certain ranges of parameters values, we suggest that a bifurcation diagram, 
the variation of the Lyapunov exponent, the existence of a periodic orbit of 
period not equal to a power of two and positive topological entropy are some 
techniques that can give a clear answer to this problem. 



1 See, e.g., [20], [18], [3], [9] and [5]. The potential for very complex behavior signif- 
icantly increases if models become somewhat less reductionist, e.g., if heterogenous 
agents and different learning processes are also taken into account, [10], [6]. 
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1 The Matching Model 



Why would the labor market in most of the developed economies behave in 
such a volatile way as the evidence points out? One possibility, and usually 
the most favoured one in the dominant view of economics, is that the economy 
has an inherently linear structure and is hit by permanent exogenous shocks. 
As these shocks are entirely unpredictable, they render the dynamics and the 
cycles hardly predictable and controllable. Another more recent view, which 
should also be considered for discussion because it seems consistent and re- 
alistic, is based on the possibility that the economy has a structure that is 
nonlinear and the cycles are an endogenous manifestation of this character- 
istic, either with or without external shocks added to the structure. In what 
follows, a very simple and fully deterministic model will be presented that is 
capable of generating such type of volatility with standard parameter sets and 
no noise. 

Let us assume that in every period of time there are large flows of workers 
moving into and out of employment: a certain number of job vacancies is 
posted by firms (v t ) and there is a total measure of workers looking for jobs 
(ut). When a worker and a firm reach an agreement there is a successful match, 
and the total number of these matches is given by the aggregate matching 
function 



The intuition behind (1) is very simple: the higher is u, the easier it will be 
for firms to get a worker with the desired qualifications; and the higher is the 
level of vacancies posted by firms v, the higher is the probability that a worker 
will find an appropriate job. For simplicity we will assume A as a constant. 
However, a more adequate treatment would consist of treating A as a variable 
dependent on the level of public provision of information by public agencies 
with the objective of increasing the number of successful matches. 

The measure of labor tightness is given by the ratio 9 t = v t /u t . Then, the 
probability of a vacancy being filled at t is given by 



Let n t+ i be the total number of employed workers at the beginning of t+ 1 and 
let s be defined as the probability of a match being dissolved at t. Therefore 
we have 

n t+l = (l-s)nt + q (9 t ) v t , 

where 9 t = v t /u t = v t /(l — n t ). Notice that (1 — s) n t gives the number of 
undissolved matches prevailing at t and passed on to t + 1, while q{9 t )v t 



M (u t , v t ) = Au1v\- a A > 0, a e (0, 1) . 
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represents the number of new matches formed at t with the available number 
of unemployed workers and vacancies. 



As shown in [12], the model can be solved for the decentralized outcome of a 
Nash bargaining game between workers and firms but to keep the model as 
close as possible to the presentation in [1] and [2] we should focus upon the 
central planner solution to the matching model. The objective function of the 
central planner is given by 

U (ra, v) = <pn t + z (1 - n t ) - cv t 

where 0, z and c are parameters that represent, respectively, the productivity 
of each worker, the lost value of leisure due to labor effort, and the cost that 
firms incur per vacancy placed in the market. 2 Therefore, the planner chooses 
v t and the next period's employment level, n t+ i, by solving the following dy- 
namic optimization problem 



^max ft [4>n t + z (1 - n t ) - cv t ] 

subject to 

n t+ i = (1 - s) n t + q [ — - ) v t , 



. 1 - ra t , 

where (5 is the time discount rate and an initial condition n is given. The 
Lagrangian can be written as 

L = [ft [<M + 2 (1 - rat) - cu t ] + A* (l-s)n t + q v ~ Ut + l } ■ 

The first order conditions (FOC), for an interior solution, are given by 
dL 



^- = -pt c + \ t W(9 t )9 t + q(9 t )]=0 



= -At + ft +1 {4>-z) + A t+1 I (1 -s) + q' (9 t+1 ) e 2 t+1 



0. 



The very interesting point in [1] and [2] was the manipulation of these FOC 
to arrive at a reduced equation that can lead to chaotic dynamics. From the 
first FOC we get At = g /(g t )g t + g (e t ) an d substituting this and the corresponding 
expression for A t+ i into the second FOC we obtain 

a9? +1 -bd t+1 = d«-d, ae(0,l) (2) 

2 Notice the trade-off between vacancies and unemployment in this objective func- 
tion. The first right hand term represents the benefits to society from successful 
matches (working) , while the last two give the leisure costs and the costs associated 
with posting vacancies. 



4 



with the following parameter definitions and restrictions 



a = (5 (1 - s) e (0, 1) , 1 > b = Aa(3 > 0, d = (A/c) (1 - a) (5 {(f> - z) > 0. 

(3) 

Equation (2) gives the law of motion for the index of labor market tightness 
in the economy under the planner's solution. In other words, given an initial 
condition 9 , equation (2) completely characterizes the trajectory of 9 and the 
whole economy. So, the backward dynamics of this model can be characterized 
by the four-parameter one-dimensional family of maps g : [0, g max ] — ► [0, 5Wx] , 
where 

g (0) = (a9 a - b9 + d)° , (4) 

i 

with the parameter restrictions in (3) and max = where g max is 

implicitly defined as the lowest positive root of the equation ag^ ax — bg max +d = 
0. The first derivative of the map g can be calculated as 

g > (6) = ( a e a -be + d) 1 ^ (V- 1 -^j,9e[0, <? max ] , 

which implies that g is unimodal with a unique maximum (critical point) at 
max . In addition, g has a unique fixed point located to the right of # max if 
g (#max) > #max- The unique fixed point of g is denoted by 6* and is implicitly 
given by a9" — = 9" — d. Despite the impossibility of the computation of an 
explicit expression for 6*, the unicity of this solution is obvious by considering 
fi (9) = a9" — b9* and f 2 (9) = 9® — <i,where fx (9) is monotonically decreasing 
for 9 from # max to +oo and f 2 (9) is monotonically increasing for 9 from to 
+oo. Therefore f\ (9) = f 2 (9) has a unique solution for 9 > 9 max . 

The fixed point is an attractor in the case of backward dynamics if \g' (9*)\ < 1 
and in forward dynamics if g' ' (9*) < — 1. For g' (9*) = — 1 a period-doubling 
bifurcation occurs and the fixed point changes stability. Since it is not possible 
to obtain a closed form expression for 9*, this condition cannot be checked in 
general but can be checked for each set of parameters separately. 



2 Chaotic Dynamics in the Model 

The first main objective is to show that the matching model leads to chaotic 
motion in the backward dynamics and the second purpose is to provide a very 
rigorous method to compute topological entropy for uni(multi)-modal maps. 
Generally, the search for chaos in one-dimensional maps is mainly based on 
the application of the well-known Li-Yorke theorem, but when a period three 
fails to be encountered it is necessary to use some other criteria. In what 
follows, we use the kneading sequence (the trajectory initiated at the critical 
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point of the unimodal map) and the associated Markov transition matrix in 
order to compute the topological entropy of the map. These are tools from 
symbolic dynamics which provide simple and exact methods to conclude about 
the existence of chaos in multimodal maps. 

In [2] the authors conclude that the unique fixed point 9* of the map (4) is 
always stable in backward dynamics, for almost any choice of the parameters 
values with respect to conditions (3). We obtain a different result: while for 
most reasonable sets of parameter values the fixed point is stable, however 
there are also other sets of parameter values that allow for unstable chaotic 
motion in the model as we show next. 

Proposition 1 For a = 0.21, (5 = 0.955, A = 0.99725, 7 = (0 - z)/c = 1.31, 
s = 0.1518, we obtain chaotic motion for the nonlinear map presented in (4)- 
Moreover, this is true for almost any value of a in the interval [0.2,0.22] . 

PROOF. The first point which we would like to emphasize is that all param- 
eter values satisfy the restrictions presented in equation (3). That is 

a = P(l-s) = 0.21 G (0, 1) , 1 > b = Aap = 0.2 > 0, 
d = A(l - a) /fy = 0.9856 > 0. 

This is illustrated in Figure 1, where the unimodal map and the typical chaotic 
time series associated to the unstable fixed point 6* = 3.6793 are presented. 
The first derivative of the unimodal map is g' = —1.8554, which shows 
that the equilibrium is unstable. By varying the parameter a in the interval 
[0.2, 0.3] we obtain the classical period-doubling route to chaos which is also 
illustrated in Figure 1. So, for any value of a arising after the 2°° bifurcation 
point, we encounter chaotic dynamics. We should point out that in order to 
obtain the bifurcation diagram illustrated in Figure 1 the parameters a, b, and 
d have also to have their values changed, as they are continuously depending 
on the values of a. The parameter restrictions are always satisfied, since a = 
0.81, 0.1904 < b < 0.2875, 0.8733 < d < 0.9980. 



Moreover, in Figure 2 the bifurcation diagram is illustrated when the parame- 
ter a is varied between 0.7 and 0.9 and taking the same parameter calibration 
as in the above Proposition: a = 0.21 ,b = 0.2 and d = 0.9868. There is 
also a period-doubling route to chaos, but the map g dynamics escapes to 
infinity before the occurrence of a period three orbit. This presents one more 
reason to consider the complete Sharkowski relation between the periodic or- 
bits of the map and to apply some other techniques in the search for chaos in 
this specific model. In the same figure we also represent the variation of the 
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Fig. 1. The map g, a typical irregular (and chaotic) orbit and the bifurcaion diagram 
when a is varied 



Fig. 2. Bifurcation diagram when the parameter a is varied and corresponding Lya- 
punov exponent. 

Lyapunov exponent for a belonging to the interval [0.7,0.9]. We recall that 
positive Lyapunov exponent is a necessary condition for chaos. 

Bhattacharya and Bunzel [1] suggest three examples for the study of the dy- 
namics of the map g. In the first case a period 3-cycle is found, which implies 
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the existence of chaos in the Li-Yorke sense if the Sharkovsky order is applied. 
In the second example, it was argued that a period 3-orbit could not be found 
for the following parameter values: a = 0.75, b = 0.58, d = 0.62, a — 0.15 since 
the equation g 3 (9) = 9 has no solution. Since the map g is unimodal and Mi- 
tra's sufficient condition for chaos in unimodal maps is verified, the conclusion 
was that for this parameter setting chaotic motion in the backward dynamics 
could also be found. 3 Finally, in the third example the set of parameters are: 
a = 0.9, b = 0.7, d = 0.6, a = 0.2. For these values no period three orbit is 
found, neither the sufficient condition of Mitra is verified. Therefore, it was 
argued that for this case the very existence of chaos for the unimodal map is 
questioned on the grounds of a lack of logical proof of such dynamics. 

In order to clarify some of the issues raised by [1], a symbolic dynamics ap- 
proach is developed for the unimodal map g, which allow us to perform the 
computation of the topological entropy for any choice of parameters, and, of 
course, permit us to classify the complexity of the map since positive topolog- 
ical entropy implies the existence of chaotic dynamics. We will concentrate on 
their third example. 

We consider again the unimodal map g : [0, g max ] — ► [0, <?max] • This kind of map 
has symbolic dynamics relative to a topological Markov partition generated 
by the orbit of the critical point 6> max . This is illustrated in Figure 3 for the 
parameter values presented in Example 2. So, any numerical trajectory 9 = 
9 niax ,9 1 ,9 2 , ... for the map g corresponds to a symbolic sequence o§0\Oi... 
where Oi G {L, C, R} depending on where the point 9i falls in, i.e., 



<Ti (00) 



L if «f (0 O ) <0* 
Cif g* (0„) = D 
RXg* (9 )>9 1: 



All symbolic sequences made of these letters may be ordered by the natural 
lexicographical order L < C < R and 

Defining the fullshift X2 = {er = t7ocricr 2 .... where cr, = L or R} to be the set 
of all possible infinite symbolic strings of L's and i?'s, then any given infinite 
symbolic sequence is a singleton in the fullshift space. The Bernoulli shift map 
s : E 2 — > S 2 is defined by s (a) = s ((7o<7i<72...) = <7i<72<73.... In general, not 
all symbolic sequences correspond to the trajectory of an initial condition 9q- 
Restricting the shift map to a subset of £2 consisting of all the itineraries that 
are realizable yields the subshift ScE 2 . 



3 As argued in [16], for a continuous unimodal map / : X — ► X, where X is a 
non-negative interval, x max is the critical point such that / (x max ) > x max and x* 
is the unique fixed point of the map such that > x max , Mitra states that: If / 
satisfies f 2 (x max ) < x max and f 3 (x max ) < x*, then (X,f) shows topological chaos. 
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Fig. 3. Markov partition for a period 5 orbit. 

We formulate the result in terms of topological Markov chains, a special class of 
subshifts of finite type where the transition in the symbol sequence is specified 
by a — 1 matrix. Any (n x n) binary matrix M = {Mij) i - =Q n _ 1 , M^- e 
{0, 1} generates a special subshift 

S M = {a e S 2 : M a . a . +1 = 1, Vi G N} 

which is called the topological Markov chain associated with the Markov ma- 
trix M. We say that M ai(Ti+1 = 1 if the transition from <jj to a i+ i is possible. 
The matrix M gives a complete description of the dynamics of the unimodal 
map. The premier numerical invariant of a dynamic system is its topological 
entropy defined by 

n->oo n 

where W n (E M ) is the set of words of length n occurring in sequences of S M . 
Moreover, if the topological Markov matrix M is given then, the topological 
entropy is the natural logarithm of the spectral radius of M. 

For the parameter values a = 0.75, b = 0.58, d = 0.62, a = 0.15, we found a 
period 5 orbit: {1.8549,0.0013,0.4756, 1.1047,0.1350} which is shown in Fig- 
ure 3 with the corresponding 4 interval {Lj} i=l 4 Markov partition. The 
critical point assumes the value # max = 0.1452 and generates the symbolic 
partition for the map g. The periodic orbit has the following symbolic ad- 
dress: {9i9 2 9- i 9 A 9 5 )' X) = (RLRRC) 00 and in consequence we have the following 
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Markov matrix: 











1 



1 















1 



Mrlrrc 







1 



1 







1 















The maximal eigenvalue of M is given by A = 1.5128, which implies that the 
topological entropy is positive: h top ~ 0.4140 and this shows very clearly that 
we are dealing with chaotic motion in this set of parameter values. 

It should be noted that for any other kneading sequence, for a suitable choice 
of parameters values, we can obtain a Markov partition and a Markov ma- 
trix which totally determine the complexity of the unimodal map. This is a 
very simple and rigorous way to estimate the topological entropy of a one- 
dimensional model and to check for the existence of chaos. 



3 Concluding Remarks 

The general use of the Li-Yorke theorem or the Mitra condition may lead 
to misleading conclusions about the existence of chaotic dynamics, as it was 
done in [1]. Therefore, in order to obtain relevant answers to whether there 
are or not chaotic dynamics under certain ranges of parameters values in a 
1-dimensional particular model, we suggest that a bifurcation diagram, the 
variation of the Lyapunov exponent, the existence of a periodic point of pe- 
riod not equal to a power of two, and symbolic dynamics are very powerful 
techniques for that purpose. Moreover, the application of these techniques to 
the matching labor market model so interestingly developed by Bhattacharya 
and Bunzel clearly confirmed that a very simple model of the labor market, 
with well behaved aggregate functions (continuous, twice differentiable and 
linearly homogeneous) do really produce chaotic behavior for a range of pa- 
rameter sets which had been questioned in [2] , for high values of the measure 
of labor market tightness. 

Financial support from the Fundagao Ciencia e Tecnologia, Lisbon, is grateful 
acknowledged, under the contract No POCTI/ ECO /48628/ 2002, partially 
funded by the European Regional Development Fund (ERDF). 
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